< 
(N 



C/3 



o 

(N 
> 

O 

o 

On 

-i— > 
c3 



i 

C 

O 

o 



X 



Normalization factors for magnetic relaxation of small particle systems in non-zero 
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We critically discuss relaxation experiments in magnetic 
systems that can be characterized in terms of an energy bar- 
rier distribution, showing that proper normalization of the 
relaxation data is needed whenever curves corresponding to 
different temperatures are to be compared. We show how 
these normalization factors can be obtained from experimen- 
tal data by using the Tln(t/ro) scaling method without mak- 
ing any assumptions about the nature of the energy barrier 
distribution. The validity of the procedure is tested using a 
ferrofluid of FeaCU particles. 

PACS numbers: 75.50.Tt, 75.60.Lr. 



I. INTRODUCTION 

The study of the relaxation of magnetic systems pro- 
vides a way to obtain information about different proper- 
ties that cannot be so easily achieved by other methods. 
Most of the works in this field are based on the logarith- 
mic or critical volume approximation!!!. To characterize 
the time dependence of the magnetization they analyze 
the relaxation rate, also called magnetic viscosity, S, as 
a function of the external parameters. When plotted as 
a function of the magnetic field it is possible to [Study 
the variation of the energy barriers with the fielda, the 
interaction effects among the .aaagnetic entities&u, and 
the switching field distributioroQ among other magnetic 
properties. When plotted as a function of the tempera- 
ture for a given field it gives information about the rever- 
sal mechanisms in films and small particle systemso and 
it has been proposed as a way to test the possibility of 
observing macroscopic quantum tunneling (MQT) effects 
at low enough temperatureaj. It is obvious that in this 
case the initial and final states of the relaxing magneti- 
zation are in general temperature dependent. Therefore, 
this dependence must be corrected in S data in order 
to compare values obtained at different temperatures. If 
this correction is not taken into account this trivial con- 
tribution to the thermal dependence of S can obliterate 
the real S(T) behaviour arising from the relevant phys- 
ical mechanism and it can even give rise to misleading 
interpretations. 

x, in some cases, Tln(i/ro) scaling has been 



usccoErl to confirm the existence of MQT by adducing 
that the fail of scaling of low temperature relaxation 
curves could be the signature of non-thermal mecha- 



nisms. As we will show later this lack of scaling could 
be only a consequence of a non-proper normalization of 
the data. 

In most of particle systems in the blocked regime and 
due to the long-time decay towards the equilibrium state, 
it is very difficult to obtain a direct experimental determi- 
nation of the equilibrium magnetization when a magnetic 
field is applied (note that, in general, the field-cooled 
state is not a true equilibrium state) . It is the purpose of 
this article to give a detailed account of the normalization 
procedure of the relaxation, curves based on the so-called 
Tln(i/ro) scaling methodOEl from which an indirect de- 
termination of the thermal dependence of the equilibrium 
magnetization (final state) can be obtained. We also dis- 
cuss the consequences of this normalization procedure 
on the analysis of the S(T) behaviour. The proposed 
method is illustrated by using experimental data from a 
ferrofluid composed of Fe3C>4 small particles. 



II. NORMALIZATION FACTORS FOR 
RELAXATION CURVES 

In the study of time dependent processes in small parti- 
cles systems two kind of experiments (which will be called 
type A and B) can be distinguished according to what 
is the final equilibrium state of the system. In Type A 
experiments the system evolves towards a demagnetized 
state in zero applied field after a previous cooling in the 
presence of a field H (FC process) and the variation of 
the thermoremanent magnetization (TRM) is measured. 
If t is the time elapsed after the field was reduced to zero 
then, in the critical volume approximations: 

/■OO 

M TRM (T, H,t)= / dEM FC (E)f(E) 

JE c (T,t) 

poo 

~M FC (T,H) / dEf(E). (1) 

JE c (T,t) 

where E c is the critical energy barrier E c (T,t) — 
Tln(t/To) which indicates the onset of superparamag- 
netic (SP) behaviour. 

In Type B experiments a zero-field cooled (ZFC) sam- 
ple increases its magnetization in a magnetic field H and 
the variation of the isothermal remanent magnetization 
(IRM) is measured. In the critical volume approximation 
the time, dependence of the magnetization is given in this 
case byfl: 
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p E ( T t H ) 

M IRM (T,H,t)= I dEM eq {E)f{E) 
Jo 

— M eq (T, H) ( dEf(E). (2) 

Jo 

In order to compare relaxation curves measured at dif- 
ferent temperatures it is necessary to remove the thermal 
dependence of the initial and final states of the magne- 
tization. It is clear from Eq. [l] and || that, due to the 
fact that in both expressions the integrals are bounded 
between and 1, this can be simply achieved by divid- 
ing magnetization data by a normalization factor which 
in Type A experiments is Mpc (T, H) and in Type B is 
M eq (T,H). 

In Type A experiments, the normalization factor, 
Mpc(T), in ([!]) comes from the contribution of the 
blocked particles to the initial magnetization or, in other 
words, from the irreversible component of the FC mag- 
netization at the temperature T. In many small particle 
systems at low T this quantity can be considered as -a 
constant in the range of temperatures usually studiedB. 
In fact, when one represents the relaxation data as a 
function of the scaling variable Tln(t/ro) all the curves 
recorded at different temperature usually superimpose 
onto a-.unique master curve without any normalization 
factoilil 

In contrast, in Type B experiments the magnetic field 
is not zero and SP particles have a temperature depen- 
dent contribution to the magnetization while blocked 
particles are randomly oriented giving no net contribu- 
tion to the magnetization in the field direction. Now the 
reversible component of the ZFC magnetization gives the 
main contribution to the normalization factor, M eq (T). 
In many cases, at low fields, M eq (T) follows a Curie- Weiss 
law and cannot be considered as a constant. The need 
of this temperature dependent normalization factors is 
clearly manifested as a vertical shift of the curves when 
data are represented in a Tln(t/ro) scaling plot. 



III. EXPERIMENTAL RESULTS 

The studied sample is a ferrofluid composed of Fe3C>4 
small particles with volume packing fraction e = 0.01 
which was obtained from chemical deposition of Fe +2 and 
Fe +3 sulfides and dispersed in an hydrocarbon oil. The 
sample analyzed by transmission electron microscopy 
(TEM) was prepared by wetting a carbon film mounted 
on an Au grid with the ferrofluid and subsequently drying 
it with air. 

In Fig. [j] an electron micrograph of the magnetic par- 
ticles is shown. Taking a sampling of 200 particles and 
considering them spherical, the particle size distribution 
have been fitted to a logarithmic-normal function with 
a = 0.24 and mean particle diameter of about 4 nm (see 
inset of Fig. II]). 



A. Magnetic characterization 

The magnetic study was performed with a commercial 
SQUID magnetometer. The highest applied field was 55 
kOe and the lowest temperature was 1.8 K. 

Figure |^ shows the thermal dependence of the magne- 
tization in an applied magnetic field of 10 Oe following 
a ZFC-FC process. The magnetic behaviour displayed is 
typical of SP particles. Above Ti rr = 50 K the system 
is in the SP region, in which the magnetization curve is 
reversible. The ZFC curve developes a maximum at Tg 
which is about 14 K. 

To obtain the mean magnetic moment of the particles 
as a function of the temperature, p-(T), we have fitted 
the magnetization curves deep in the SP region (within 
the range 80 to 200 K) to a distribution of Langevin 
functions following the procedure described in Ref. Kq. 
By extrapolating these values to T = K with a T^i 2 
law (as expected according to a spin-wave theory for a 
ferrimagnetic materials), a mean magnetic moment per 
particle of /2 = (1320 ± 20)/is is obtained. Using j2 and 
the bulk magnetic moment of the magnetite the mean 
magnetic size of the particles has been found to be 3.5 
nm in diameter. 

Tb has been measured as a function of the field. By 
extrapolating these values to Tb = we have estimated 
the mean anisotropy field (the field at which the mean 
energy barrier disappears) to be Ha — 5000 Oe. From 
Ha the mean anisotropy constant of the particles have 
been estimated as K — JIHa/^V = 1.3 x 10 5 J/m 3 , where 
V is the mean particle volume. This value is fik>f~i to 
others found in the literature for similar systemsOJilJ. 

For the subsequent analysis of the relaxation curves 
it is also convenient to study the variation of the ini- 
tial susceptibility with the temperature. For a system 
of interacting magnetic particles in the SP regime this 
quantity can be written astj: 



3k B (T-T ) 



(3) 



where To is due to the existence of interaarticle interac- 
tions and/or the effect of f(E) (see Refill). The inset in 
Fig. H shows the inverse susceptibility as a function of the 
temperature as obtained from the ZFC-FC curve, where 
the thermal variation of the mean moment of the parti- 
cles has been corrected by using the T 3 / 2 law obtained 
before. Above Ti rr ,l/x shows a linear dependence on 
T. Fitting 1/x to (g) we have estimated T = -11 ± 3 
K. This value of To may be mainly attributed to the 
existence of demagnetizing interparticle interactions tak- 
ing into account that the volume distribution is not very 
broad. 
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B. Normalization factors and Tln(tro) scaling 

The magnetic relaxation curves at different tempera- 
tures were recorded with a SQUID magnetometer follow- 
ing the procedure described in Ref. |l9|. The measure- 
ments were started 5 s after applying the field and were 
performed during aprox. 1000 s at the lowest tempera- 
ture and aprox. 10000 s at the highest temperature. The 
relaxation curves after ZFC the system were measured in 
the presence of a magnetic field of 10 Oe (type B exper- 
iment) while in those measured in zero field the system 
was previously FC in 10 Oe (type A experiment). 

In the following we will describe an experimental pro- 
cedure, based on the Thx(t/ro) scaling methodlI3, to ob- 
tain normalization factors, M eq (T), for relaxation data 
recorded in the presence of a field since they cannot be 
directly measured due to the long-time decay of the mag- 
netization. 

Firstly, the attempt time To = 3 x 10~ n s has been 
evaluated by scaling the relaxation curves in zero field 
after FC the system at several temperatures following 
the method previously described in Ref. 12. For this 



purpose we have considered that the initial magnetiza- 
tion Mpc{T, H) is constant in the temperature range of 
the experiment so no normalization factors are needed to 
scale the curves. We will assume that the field variation 
of t q is smaller than the error in itSpdetermination at low 
fields, according to Brown's theoryc2l. 

In the next step, relaxation data recorded in the pres- 
ence of a field have been plotted as a function of the 
scaling variable using the value of To previously deduced 
(see Fig. |^). After this transformation the relaxation 
curves at different temperatures are separated along the 
vertical axis by temperature dependent shifts. Taking 
into account that the applied field is much lower than 
Ha and therefore the energy barriers have not been sig- 
nificatively affected, this lack of scaling is a clear demon- 
stration that magnetization data must be normalized to 
achieve an equivalent scaling to that obtained in zero 
field. As has been discussed in Sec. [n] the normaliza- 
tion factors are proportional to M eq (T, H). If now we 
assume that M eq (T,H) are-given by Mfc(T, 10 Oe), as 
suggested by some authorsll3, no scaling is achieved be- 
cause FC magnetization does not correspond to the true 
equilibrium state. Note that M eq (T,H) can not be cal- 
culated without making any a priori hypothesis about 
the form of f(E) and the magnetic microstructure of the 
system. 

The normalization factors can be found by referring 
the different curves to the lowest temperature one. Once 
this process has been performed, the relaxation curves 
collapse onto a unique master curve that is shown in Fig. 
|]. The values of the normalization factors follow a Curie- 
Weiss law of the form (||) with To = — 15±2 K (see inset of 
Fig. |^). The extrapolation of this law superimposes with 
the susceptibility corresponding to the reversible (SP) 
region obtained from ZFC-FC measurements (see Fig. ||), 



reflecting the fact that for long enough observation times 
all the particles have become SP and demonstrating that 
the normalization factors are proportional to M eq (T, H). 



Magnetic viscosity and energy barrier 
distribution 



The classical magnetic viscosity commonly defined as 
S(t,T) — dM(t)/d(\nt), cannot be directly compared 
at different temperatures because: 1) usually magneti- 
zation is not normalized (initial and final states of the 
relaxation process change as the temperature varies) and 
2) if relaxation data have been recorded in a fixed time 
window the energy barriers which are relaxing at dif- 
ferent temperatures are not the same. Both problems 
can be circumvented by defining magnetic viscosity as 
S(t,T) = dM(t)/d(T\n(t/T )) (M(t) is the normalized 
magnetization used in the scaling procedure). By per- 
forming the Tln(t/ro) derivative S measures the relax- 
ation rate of the magnetization due only to the energy 
barriers around T\n(t/rp). On the other hand, as pre- 
viously noted in Ref. ElJ, S is a magnitude proportional 
to the energy barrier distribution and therefore it has a 
direct physical meaning. 

These two magnitudes are simply related by 



S 



M FC T 

in the case of type A relaxation experiments, and by 

S 



s = 



M eq T 



(4) 



(5) 



for type B experiments. In the first case, and for sys- 
tems with a certain degree of interaction between parti- 
cles Mpc is usually almost temperature independent and 
both magnitudes differ by a T _1 factor. Therefore con- 
clusions from the thermal variation of S obtained from 
type A experiments must be carefully derived. Note in 
particular that if S happens to be temperature indepen- 
dent in a certain range, a result which could be inter- 
preted as aj-nroof of the existence of quantum relaxation 
phenomenal^, this would be a consequence of an energy 
barrier distribution proportional to 1/E in this range, as 
the S ~ thermal dependence reveals. It is worth 
noticing that if there is not a certain degree of freezing 
due to magnetic interactions, Mpc can not be considered 
as a constant and its thermal variation must be corrected 
in S. However, in type B experiments, only if M eq (T) is 
inversely proportional to the temperature (Curie law), 
as is the case for a sample with no or very small inter- 
particlediiteractions and a narrow distribution of energy 
barriersE£l, both magnitudes nearly coincide (note that 
this is not the case of the sample studied in this paper) 
because the thermal variation of M eq cancels the factor 
T in Eq. g. 
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For the sample studied in this paper 5(Tln(f/ro)) has 
been obtained by performing the numerical derivative of 
the master curve of Fig. ^ and has been compared to the 
viscosity S(T) as obtained from the logarithmic deriva- 
tive of the relaxation data at each temperature. The 
results are shown in Fig. || together with the energy 
distribution obtained from the electron micrography by 
expressing the volume distribution in energy units with 
the help of the value of K previously derived. The co- 
incidence between S and f(E) shows the consistence of 
the normalization used in the scaling procedure for type 
B experiments. On the other hand, it is important to 
note that S does not coincide with S because for this 
sample M eq is not simply proportional to T . Only the 
overall shape of the energy barrier distribution obtained 
from TEM is reproduced by S, but shifted to higher ener- 
gies. Note also that the extrapolation of the quasi-linear 
low temperature regime of S intercepts the temperature 
axis at a non-zero value as has been reported in other 
systems. In our case, this result is only a consequence of 
the lack of normalization and has no physical meaning. 

IV. CONCLUSIONS 



low T if MQT occurs. If this is the case, even a multi- 
plicative factor is not enough to superimpose relaxation 
curves obtained at different T (T is then an irrelevant pa- 
rameter that would not have to be included in the scaling 
variable). 

In conclusion, we have shown that the Tln(t/ro) scal- 
ing method provides a useful tool to obtain the normal- 
ization factors and the energy barrier distribution in both 
kinds of experiments and even in systems that cannot be 
considered as an assembly of independent small particles 
(i.e. multilayered systems, cluster glasses, amorphous al- 
loys, etc.) without making any assumptions about the 
nature of the magnetic microstructure. 

ACKNOWLEDGMENTS 

We are indebted to professor S.W. Charles for provid- 
ing the sample studied in this paper. Finantial support 
from both the Spanish CICYT through MAT94-1024- 
CO2-02 and the Catalan CIRIT through GRQ1012 is 
acknowledged. 



We have stressed the importance of proper normal- 
ization whenever relaxation curves measured at different 
temperatures must be compared. In the case of exper- 
iments performed in zero field, care must be taken in 
systems for which Mpc(T) cannot be considered as a 
constant. When this is the case, non-normalization could 
give place to a spurious thermal dependence. 

In the case of relaxation experiments performed in an 
applied magnetic field, there exists a certain controversy 
in the literature about the nature of the normalisation 



: either no normalization factors are used at alflc3l2J or 
the FC magnetization value, corresponding to the field 
at which the experiment is performed^, is used. 

In systems for which the M eq (T) follows a Curie law 
(non-interacting particles, neglictible To) the first option 
happens to be correct by chance as can be easily seen in 
(||). In this kind of systems the second option is partic- 
ularly wrong when applied to low temperature data be- 
cause the FC magnetization is slightly temperature de- 
pendent while the SP magnetization, which is the true 
equilibrium state at long times, follows a T _1 behaviour. 

In fact, when Tm(f/ro) scaling is used to evidence 
quantum relaxation mechanisms through a lack of scaling 
of relaxation data (recorded in the presence of a field), 
no conclusions should be extracted without previously 
having tried to norm alize data following the process de- 
scribed in Sec. 



[II B 



The sample studied in this paper 
is an example where a clear fail of scaling of the non- 
normalized data does not indicate any non-thermal pro- 
cess (S agrees with the energy distribution deduced from 
f(V), see Fig. ^|). This does not mean that it is always 
possible to find scaling factors for data corresponding to 
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FIG. 5. Reversible region of the reciprocal susceptibility 
(open squares) and thermal dependence of the inverse of the 
normalization factors (solid circles) necessary to join the re- 
laxation data of Fig. ^ onto a unique master curve. Solid 
line is a linear regression of both data. Inset: inverse of the 
normalization factors as a function of the temperature. Nor- 
malization factors have been reduced to the same units of the 
susceptibility multiplying them by an arbitrary quantity. 

FIG. 6. Effective distribution of energy barriers as obtained 
from the numerical derivative of the master relaxation curve 
(solid line). Magnetic viscosity S(T) as obtained from the 
logarithmic time derivative of the relaxation data at the tem- 
peratures indicated in Fig. ^ (solid squares). The energy 
distribution f(E) obtained from f(V) is also shown for com- 
parison (open circles). 



FIG. 1. Electron micrography of the sample obtained by 
TEM. The inset shows the distribution of particle diameters 
obtained from a sampling of 200 particles. The solid line is a 
logarithmic-normal function with a — 0.24 and mean particle 
diameter of about 4 nm. 



FIG. 2. Temperature dependence of the magnetization of 
the sample for a ZFC (lower curve)- FC (upper curve) pro- 
cess in a magnetic field of 10 Oe. Inset: reciprocal of the 
ZFC-FC susceptibility corrected from the thermal variation 
of the mean particle moment. 



FIG. 3. Relaxation data recorded from 1.8 to 15 K in the 
presence of a magnetic field of 10 Oe after ZFC the sam- 
ple as a function of the Tln(t/ro) scaling variable. Open and 
solid simbols correspond alternatively to the temperatures in- 
dicated in the Figure. Inset: detail of the lowest temperature 
region. 



FIG. 4. Scaling plot for the relaxation measurements 
shown in Fig. |j| Open and solid simbols correspond alter- 
natively to the temperatures indicated in the figure. Inset: 
detail of the lowest temperature region. 
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